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ABSTRACT 

We use the operator product expansion (OPE) to show that non- 
perturbative QCD corrections to Ap can be calculated using unsub- 
tracted dispersion relations for either the transverse or the longitu- 
dinal vacuum polarization functions. Recent calculations of the non- 
perturbative contribution to Ap based on a non-relativistic calculation 
of corrections to the tt threshold are inconsistent with this result. 
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I. INTRODUCTION 



An accurate calculation of the contribution of the t—b doublet to the p pa- 
rameter is important in two respects. The first is that it will let us translate the 
tight constraint placed on the p parameter by LEP measurements to limits on 
the t quark mass, and the second is that in the event that the t quark mass is 
measured directly at the TEVATRON, it will help us constrain the contribution 
of new physics to Ap. 

The contribution of the t—b doublet to the p parameter has been calculated 
to 0(aa s ), and the result is given by 



Ap 



3am 2 



16irs 2 c 2 m 2 z 




(i; 



in the limit — > [|T], 0. 

Recently, Kniehl and Sirlin estimated the size of the higher order QCD correc- 
tions to Eq. (JJ) using a dispersion relation for Ap [[J. Their approach has been 
to assume that the effect of non-perturbative QCD corrections is dominated by 
the change in the shape of the ti threshold. This change is then calculated using 
the leading non-relativistic approximation, and substituted into the dispersion 
relation for Ap. This work has lead to a certain amount of controversy since 
different authors found the effect to be different in magnitude, ranging from 10% 
to 80% of the 0(aa s ) correction, and sometimes even different in sign ||. ||. 

The difference in sign comes from the fact the there are two possible dispersion 
relations for Ap. In order to explain what they are, we must introduce some 
notation. Following Ref. [0], we define 



Uli A (q, mi ,m 2 ) = -ij d 4 xe^(0\T* [jJ' A (x)JJ' A \0)\ |0> 
= g fil/ U v ' A (s, mi, m 2 ) + q^q v \ v ' A {s, m 1 , m 2 ) 



n^(s,mi,m 2 ) + — — A^(s,mi,m 2 ; 



9fJ,f > 

v T 



(2) 



where s = q 2 , and JY' A (x) represents the vector and axial vector currents con- 
structed from the fermion fields, respectively. Note that 



so that 



A v ' A (s) = U v ' A (s) + s\ v ' A (s) 
A V ' A (0) = 11^(0), 



(3) 
(4) 



2 



unless A ' (s) has a pole at s = 0. We further introduce the notation 



an* 

X V ± A (s 



Tl v ' A (s,m l} m 2 ) 
\ v ' A (s,m 1 ,m 2 ), 
A v ' (s,m l ,m 2 ) 



U v,A (s, mi, mi) + U v ' A (s, m 2 , m 2 



V,Ai 



A ' (s,mi,TOi) + A ' (s,m 2 ,m 2 ) 
A v ' A (s, mi, mi) + A v ' A (s, m 2 , m 2 



The conservation of the neutral vector currents implies the Ward Identities: 



-s\ (s 



K(s) = 0. 



(5) 



(6) 



These definitions let us write the contribution of an SU(2) fermion doublet, with 
masses mi and m 2 , to the vacuum polarizations of the W and the Z at zero 
momentum transfer s = as 



rw(o) 
iWo) 



r 



n£(o) + n£(o) 



r 



g* + g 



itf(o) + itf(o) 



A^(0) + A^(0) 
A^(0). 



(7) 



Note that 11^(0) is actually zero from current conservation (c/. Eq. (||)) but we 
will keep it in our expressions for later convenience. Inserting Eq. ([!]) into the 
definition of the p parameter, we find 



Ap 



iW(o) n zz (o) 



M 2 W 



= -| { [nl(o) + n£(o)] - [rtf(o) + itf (o)] } 

= ^{[A^(0) + A^(0)]-A^(0)} 
Applying the 'unsubtracted' dispersion relation: 



f(s) = - Tds^ 

7T J S 



, Imf(s') 



s + ie 



(9) 
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to the expression for Ap using the IT(s)'s gives us 



Ap = J°° ^ [{imn^(s) + Imn^(s)} - {imn^(s) + Imll^(s)}] . (10) 

This dispersion relation is equivalent to that introduced in Ref. ||, and used in 
Ref. @ to calculate the 0(aa s ) correction. 

If we write down an unsubtracted dispersion relation for the expression of Ap 
using the A(s)'s, we get 



Ap =^lf^ [{imAl(s) + ImA^( S )} - ImA^s) 



TP 



This is the relation that was introduced in Ref. ||. Eqs. ( |10"D and ([TT|) are the 
two dispersion relations for Ap that have appeared in the literature. Our point 
here is to illustrate that, although Eqs. (|10| , |TTf ) should yield the same estimate 
for Ap, under certain commonly used approximations they do not. 

In Sec. II we clarify the conditions that the IT(s)'s and A(s)'s must satisfy for 
Eqs. fllOD and (|TT|) to be true. These conditions are indeed satisfied to order aot s 
in perturbation theory. In Sec. Ill, we argue that the operator product expan- 
sion (OPE) of the current-current correlators suggests that these conditions are 
satisfied for non-perturbative QCD corrections so that both Eqs. ([10]) and (|TTD 
are valid. In Sec. IV we show that the two dispersion relations Eqs. (0) and 
([TT|) give different answers when used in a threshold approximation to calculate 
non-perturbative effects. We argue that the disagreement between the two dis- 
persion relations can be understood as a result of neglecting the non-threshold 
contribution of the Imil(s)'s and ImA(s)'s to Ap. Furthermore, the magnitude of 
the discrepancy between Eqs. (|T0| , |TTD can be used as a measure of the accuracy 
of the calculation. We conclude in Sec. V. 

II. THE DISPERSION RELATIONS FOR Ap 

Let us begin this section by recalling that the analyticity of a vacuum polar- 
ization function f(s) and Cauchy's theorem tell us that 

m ,if d ,j^L + j_ f u m. (12) 

Ti J s' — s + ie 2m J\ S \=A 2 s' — s 

Hence for Eq. @ to be true, f(s) must satisfy 

li m I jQfl = o. (13) 

A 2 ^oo/| s |=A2 s' - S 
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Applying Eq. (0) to the 11(0) 's in Eq. (§), we find 



n£(o) 

n£(o) 
n y (o) 
itf(o) 



7T 



A 2 



Imn^(s) + — i — lilts), 
s ±K J 2ni J\s\=a* s ±K h 



ds. 



1 /"A 2 ds 



ds. 



s ±K J 2ni J\ s \=a 2 s ±K J 

1 r A2 ds T rrV , . 1 /■ ds,,, . 

- / — ImnJf(s) + * — njf (s), 

- / A2 ^i m n^( s ). + — f -n A ( s ). 



(14) 



We regard the IPs in these relations to be regularized Q and thus finite quantities 
so that both the left and right hand sides of these equations are well defined. 
Note that the Ward Identity, Eq. (0), ensures that 



n y (o) 



s 

rA 2 



2m /|s|=a 2 s 

1 r AA 1 r 

-J dsImX (s) + — f^ A2 ds\ 



0. 



(15) 



as required. This result follows simply from Cauchy's theorem and the (assumed) 
analyticity of A(s). If we define 



Ap T (A 2 ) 
AR T {A 2 ) 



G F l r A2 ds 
G p 



J -j- [{imn^(s) + Imn^(s)} - {imn^(s) + Imn^(s)} 



I ^[{Hl(s)+Ui(s)}-{U^S)+Ut(s)} 

J \s\=A S L J J 



y/2~2ni J\s\=a 2 s 



the substitution of Eq. ([L4]) into Eq. (H) gives 

Ap = Ap T (A 2 ) + AR T (A 2 ). 



(16) 



(17) 



(The subscript "T" stands for "transverse".) Therefore, for the dispersion relation 
Eq. ( |T0"D to be valid, we must have 

lim AR T (A 2 ) = 0. (18) 

A 2 ^oo 



With MS or some other regularization scheme which respects the symmetry of the theory 
so that the Ward Identities are satisfied. 
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This requires the linear combination of the n(s)'s in the integrand of Ai?^(A 2 ) 
to vanish as \s\ — > oo. 

Similarly, applying Eq. ( |T2"| ) to the A(s)'s gives us 



Al(0) 
A£(0) 



l r A2 ds 

71 J S 
1 /"A 2 rfs 
7T J S 

rA 2 d s 



1 / (is 

27T2 J\s\=A 2 S 

If ds 
lixi J\s\=A 2 s 

ds 



-ImAl(s) + i ^-Al(s), 



—ImAi(s) + ^ <j> —Ai{s), 



(19) 



Again, we consider the A(s)'s to be regularized and finite quantities. If we define 

Ap L (A 2 ) ee ^|i/ A2 ^[{lmAl( S ) + ImA^( S )}-ImA^( S ) 
Gf If ds 



Afii(A ! 



y/2 2ni 7|s|=a 2 s 



the substitution of Eq. (|19|) into Eq. (f|) gives us 

Ap = Ap L (A 2 )+Ai? L (A 2 ). 



(20) 



(21) 



(The subscript "L" stands for "longitudinal".) Therefore, for the dispersion re- 
lation Eq. (PD to be valid, we must have 



lim AR L (A 2 ) = 0, 

A 2 ^oo 



(22) 



which requires the linear combination of the A(s)'s in the integrand of Ai?x,(A 2 ) 
to vanish as \s\ — > oo. 

The actual asymptotic forms of the II(s)'s and A(s)'s up to order aa s can 
by found in Ref. |7j and one can explicitly check that Eqs. fll8D and (|22|) hold. 
Therefore, to order aa s in perturbation theory, 

Ap = Ap T (oo) = Ap L (oo), (23) 

III. THE OPERATOR PRODUCT EXPANSION 

In the previous section, we have seen that 

AR T (oo) = Ai? L (oo) = (24) 
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implies 

Ap = Ap T (oo) = Ap L (oo), (25) 
and that Eq. (|24]) holds to order 0(aa s ) in perturbation theory. Whether these 



equations continue to be valid to higher orders in perturbation theory, and non- 
perturbatively is a difficult problem in general. However, since we are only in- 
terested in QCD corrections, the non-perturbative asymptotic behavior of the 
vacuum polarization functions IT(s) and A(s) as \s\ — > oo can still be extracted 
from their operator product expansions. 

The OPE's for n± i0 (.s) and A-b >0 (s) can be found in the appendix of Ref. ||. 
They are: 



AK-Q 2 ) 
A^(-Q 2 ) 

and 
where 



C A i(Q) IMQ) ~ rn2(Q)f + C A2 (fi) [rhi(ji) - m 2 (/x)] 2 + O 



,Q 2 



Cai{Q) IMQ) + ^2{Q)f + C A2 (fx) [mi(/i) + m 2 (p)] 2 + O — 



C A1 {Q) [2m 1 (Q) 2 + 2m 2 (Q) 2 
+C A 2(fJ>) f2mi(/i) 2 + 2m 2 (/i) 2 



nS(-Q 2 ) = Q 2 A^(-g 2 ) + A^(-g 2 
[MQ) + MQ)? 



Q 2 

(26) 
(27) 



CAi(g) + CA 2 (g)- 



Q 2 



A^-g 2 



c f Ai(g) + c , A 2 (g) 



[mi(g)-m 2 (g)] s 

g 2 



+ CA 3 (g) ^ (g) ^ 2(g)]2 + o 



A y (-g 2 ) 
A^(-g 2 ) 



c , Ai(g) + c , A2 (g) 



g 2 

[2m!(g) 2 + 2m 2 (g) 2 ] 

g 2 



(2f 
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The short distance physics is embodied in the Wilson Coefficients C*(Q) which 
can be calculated perturbatively, and their explicit forms for the first few orders 
in a s (Q) can be found in Ref. ||. The long distance non-perturbative physics is 
embodied in the running masses mi i2 (Q) and the vacuum expectation values of 
the higher dimensional operators that appear at higher order in the OPE. The 
Wilson coefficients are independent of the masses except through the running of 
a s (Q). 

That only these combinations of running masses appear at dimension 2 in 
the OPE can be understood as follows: Consider the charged channel functions 
A± A (s) and A± A (s). Since they must be symmetric under the interchange mi «-> 
fn 2 , they can only depend on (mi ±ro 2 ) 2 . Changing the sign of one of the masses 
will interchange the vector and axial-vector cases so the coefficient of (mi ±m 2 ) 2 
in the vector channel is equal to the coefficient of (mi =p m 2 ) 2 in the axial- 
vector channel. Since A±(s) must vanish when mi = m 2 , it can only depend on 
(mi — m 2 ) 2 , which in turn means that A±(s) can depend only on (mi + m 2 ) 2 . The 
dependence of the neutral channel functions A^' (s) and Xq' A (s) on the running 
masses follows trivially from that of the charged channel functions. 

Note that the Wilson Coefficients C*(Q) and the running masses rhi^iQ) 
depend only logarithmically on Q. Therefore, though these OPE's are derived 
in the deep Euclidean region — s = Q 2 ^> 0, we can expect the dependence of 
the A(s)'s and A(s)'s on the powers of s to be the same all around the circle 
\s\ = A 2 . We can see immediately that this implies ARl(oo) = ARt(oo) = 0. 
We can therefore conclude that Eq. (|25|) is correct even when higher order and 
non-perturbative QCD corrections are taken into account. 

IV. CALCULATING it THRESHOLD EFFECTS 

Let us now turn to the problem of using dispersion relations to calculate the 
non-perturbative QCD corrections to Ap. In order to make use of the dispersion 
relations Ap = Ap T (oo) and/or Ap = Ap L (oo), we need to know the spectral 
functions Imll(s) and/or ImA(s) when non-perturbative corrections are taken 
into account. Since it is impossible to calculate them exactly, this means that we 
must make some assumptions and approximations about their non-perturbative 
behavior. 

Let us denote the difference between the non-perturbative and perturbative 
vacuum polarization functions, by 

Slmll(s) = ImlLvp(s) — Imllp(s), 

5ImA(s) = ImAivp(s) - ImAp(s), (29) 
In Ref. 0, it was assumed that the most important effect of non-perturbative 
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QCD corrections on ImIT(s) and ImA(s) is to modify the shape of the tt thresh- 
old. The threshold region is also the place where higher order corrections in a s 
can be resummed in the leading non-relativistic approximation and calculated 
reliably in terms of a simple non-relativistic Schrodinger Green's function, as in 
Ref. ||. Far from the threshold region, the Imil(s)'s are assumed to be well 
approximated by their 0(aa s ) perturbative results. Therefore, in this approxi- 
mation, the functions 51mll(s) and 5ImA(s) have their support only in the region 
near the threshold. Furthermore, in the leading non-relativistic approximation, 
only the s-wave states contribute so that 



^-ImIT y ( S ) = -s5 t MK(s) = -s6itlm)4(s) = -<J tt -bnA£(s) 



(30) 



while all the other <5 4 jlmn(s)'s, (%ImA(s)'s, and 5 t fImA(s)'s are zero. (We have 
added a subscript to 5 to indicate that since we don't expect large threshold cor- 
rections in the bb or tb channels, only tt threshold effects are included. This would 
not be the case if we were considering threshold effects for say a heavy fourth 
generation of quarks, but including the other channels does not alter our con- 
clusions) . Note that the first equality in Eq. (|30| ) comes from the Ward Identity 
Eq. (^|), the second equality comes from the spin independence of QCD interac- 
tions in the non-relativistic limit, and the third equality comes from Eq. (Q), and 
the fact that 5 4 flmll^(s) = in this approximation. We will not specify what 
these non-zero terms look like in any detail since it is irrelevant to the following 
discussion. 

Now, let us see what happens when we substitute Eq. ([3(]) into the definitions 
of Apr(oo) and Apl(oo). We find: 



8 t t[Ap T (oo)} = j= 



G f 

72 



1 r -feIn<( S ) 
n J s 



(31) 



and 



5u[Ap L (oo)] 



G p 

'V2 



1 r°° ds 



71 



G p 




= 75 


71 J 



r°° ds 



6ulmU^(s) 



--X, 
(32) 

which are of the same magnitude but opposite in sign. This is the disagreement 
in sign that was mentioned in the introduction. 

There are two ways to interpret this result. The first is that either one, or 
both of the dispersion relations Eq. ([[(]) and (|TT|) are wrong. This means that 
Eqs. (|l8| ) and (|22|) cannot be both correct. This is the approach adopted by the 
authors of Ref. who, for a number of reasons, prefer Eq. ( pTTj) to Eq. (|10"D. 
However, we do not think this possibility is very likely as we don't see where the 
OPE argument of the previous section could have failed. 
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A second and more plausible possibility is that using only the leading non- 
relativistic limit to calculate non-perturbative contributions to Ap is simply not 
a good approximation and that the disagreement between Eq. ([H]) and (^2j) is a 
reflection of that fact. 

Let us take a closer look at where this apparent inconsistency between the 
two approaches comes from. The reason we get the same magnitude but opposite 
signs in Eqs. (|3"ID and fl32D is because of Eq. (|30|) , which is only true at leading 



order in the non-relativistic approximation in a small region near the threshold. 
If we could increase the range where we can calculate non-perturbative effects, 
Eq. p0| ) will not be true over the entire integration interval and we expect the 
difference between the two results to be reduced. 

We illustrate this with an example: Consider the non-relativistic limit of the 
vacuum polarizations calculated to one loop in the ti channel, 



S 167T 



0(F 



ImA^,,-(s) 33 

= ^T + W > ( 33 ) 

S l07T 



where 3 — yl — Amf/s. So, if we decided to calculate the perturbative con- 
tribution to Ap using the leading non-relativistic approximation (and including 
all the other channels), we would again get answers of the same magnitude but 
opposite sign depending on whether we chose to calculate using the n(s)'s or 
the A(s)'s. However, from the full perturbative calculation, we know that both 
techniques should give the same answer. For similar reasons, (although there 
are important differences in the perturbative and non-perturbative cases), we be- 
lieve pushing the cutoff further and further away from the threshold would cause 
the non perturbative results obtained using either the II(s)'s or the A(s)'s to 
converge towards each other, though perhaps very slowly. (The importance of 
the contribution from regions away from the threshold has already been noted 
in Ref. 0.) We are currently studying how much of this convergence can be 
achieved by taking the non-relativistic approximation to higher orders in 3 and 
thus expanding its region of applicability. 

In any approximation, a good way to test its accuracy is to see how well it 
reproduces a known result. In the present case, a good way to test how well our 
approximation gives the correct value for Ap-r(oo) or Api(oo) is to see how well 
it reproduces the known value for Ap<r(oo) — Api(oo), namely zero. When seen in 
this context, Eqs. (ftH 0) are an indication that the ti threshold approximation 
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fails, indeed we cannot even determine the sign of the additional contribution to 
Ap. 

Of course, the bright side of it is that now we may understand the reason why 
the two dispersion relations Eq. (|T0|) and Eq. (|ll|) give seemingly contradicting 
results for the tt threshold effects. In fact, they are not contradicting at all. The 
difference between the two values is the error that should be associated with the 
approximation. 

V. CONCLUSION 

We have used the OPE to show that it is equally valid to calculate Ap using 
either the dispersion relations of Ref. |J, or those of Ref. 0. These correspond 
to using unsubtracted dispersion relations for either the transverse or the lon- 
gitudinal parts respectively of the vacuum polarization functions that appear in 
the expression for Ap, and includes the consideration of non-perturbative effects. 

When the two dispersion relations are used to calculate the effect of the tt 
threshold on Ap, they give results which are equal in magnitude but opposite 
in sign. This disagreement should not be interpreted as a sign that one of the 
dispersion relations is wrong, but as a sign that neglecting the non-threshold 
region when calculating non-perturbative effects to Ap is a poor approximation. 
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